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Abstract 

Background: In the environment of high neutrino-fluxes provided in core-conapse supernovae 
or neutron star mergers, neutrino-induced reactions with nuclei contribute to the nucleosynthe- 
sis processes. A number of terrestrial neutrino detectors are based on inelastic neutrino-nucleus 
scattering and modeling of the respective cross sections allow predictions of the expected detector 
reaction rates. Purpose: To provide a self-consistent microscopic description of neutrino-nucleus 
cross sections involving a large pool of Z = 8 — 82 nuclei for the implementation in models of 
nucleosynthesis and neutrino detector simulations. Methods: Self-consistent theory framework 
based on relativistic nuclear energy density functional is employed to determine the nuclear struc- 
ture of the initial state and relevant transitions to excited states induced by neutrinos. The weak 
neutrino- nucleus interaction is employed in the current-current form and a complete set of transition 
operators is taken into account. Results: We perform large-scale calculations of charged-current 
neutrino-nucleus cross sections, including those averaged over supernova neutrino fluxes, for the set 
of even-even target nuclei from oxygen toward lead {Z = 8 — 82), spanning = 8 — 182 (OPb 
pool). The model calculations include allowed and forbidden transitions up to J = 5 multipoles. 
Conclusions: The present analysis shows that the self-consistent calculations result in consid- 
erable differences in comparison to previously reported cross sections, and for a large number of 
target nuclei the cross sections are enhanced. Revision in modeling r-process nucleosynthesis based 
on a self-consistent description of neutrino-induced reactions would allow an updated insight into 
the origin of elements in the Universe and it would provide the estimate of uncertainties in the 
calculated element abundance patterns. 
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I. INTRODUCTION 



Neutrino-induced reactions on nuclei play an important role in nuclear astropyhsics, in 
particular during core-collapse supernova evolution and nucleosynthesis jl, 2|. Elastic neu- 
trino scattering on nuclei and nucleons determines the neutrino trapping and the diffusion 
time scale of the outwards streaming neutrinos jsl- A variety of processes contribute to the 
energy loss in stellar interiors, e.g. pair, photo-, plasma, nucleon-nucleon bremsstrahlung, 
and recombination neutrino processes 141. In addition, inelastic neutrino-electron scattering 

n 

mainly thermalizes the neutrino spectra [5| and inelastic neutral-current neutrino-nucleus 
scattering has a large contribution to the opacity js, 7|, during the core collapse and subse- 
quent explosion phase. Additionally, in the environment of exploding massive stars, which 
has long been explored as a possible site for the r-process nucleosynthesis, charged- current 
neutrino-nucleus reactions play an important role in the production of chemical elements. 
As pointed out in Ref. js], there are many interesting effects of neutrino-induced reactions 
before, during, and after the r-process. To explore their role in nucleosynthesis reliable 
neutrino-nucleus reaction rate compilations involving also nuclei with large neutron excess 
are required. Available data on neutrino- nucleus cross secti^s are limited to deuterium j^. 



and and ^^Fe target nucleL obtained by the LSND [10| and KARMEN [ill, [l2j col- 
laborations, and at LAMPF [l3|. Therefore, only theoretical approaches can provide cross 
sections for a large number of target nuclei that are involved in various applications of neu- 
trino physics and astrophysics. Modeling neutrino-nucleus interactions is also important 
in view of the current research and development of neutrino detectors, e.g., for supernova 
and solar neutrinos, neutrinos produced in laboratories, and geoneutrinos. The ongoing 



and planned neutrino detector facilities involve a variety of target material 



tions and scientific objectives, e.g., MOON Il4|, MiniBooNE 



OPERA 



15|, MINOS 



s, induced reac- 



161, SNO+ (l7i, 



18|, LVD (Large Volume Detector) [19|, ORLaND experiment proposal at the Spal- 



lation Neutron Source (SNS) 20|, NOvA neutrino experiment 2]J|, Daya Bay reactor neu- 



trino experiment 



22l |. etc. There is also interesting concept o 



23 



Deta-beams for the production 



24l |. The proposal to establish 



of neutrinos by using /3-decay of boosted radioactive ions 
a beta-beam facility that could produce low-energy neutrino beams in the 100 MeV energy 
range would allow direct insight into the neutrino-induced reactions in nuclei as well as its 
underlying structure involving a non-trivial combination of nuclear allowed and forbidden 
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transitions 



2M- 



At low neutrino energies, the neutrino-induced reactions are sensitive to the properties of 
nuclei involved, i.e. their initial and excited states. Therefore it is necessary to employ a mi- 
croscopic framework providing reasonable description of nuclear structure properties. Over 
the past years, a variety of advanced microscopic models have been developed and employed 
in studies of charged-current neutrino-induced reactions at low energies, also including var- 



ious particle decay channels. In particu 
random phase approximation (RPA) js]. 



ar, these include the nuclear shel 
sJaTi, continuum RPA (CRPA) 



model 



38 



25|-|311, 



421, and the 



hybrid model which combines the shell model for allowed transitions, with the RPA to ac- 
count for the forbidden transitions jsol . |43|-45|. Frameworks based on quasiparticle RPA 



(QRPA) have also recently been developed, based on Skyrme functionals 34,|46|, Brueckner 



G matrix employed for two-body interaction by solving the Bethe-Salpeter equation using 



Bonn-CD potential 



47 



49|, andprojected QRPA 50N52|. The Fermi gas model 



53 



55 



56| 



and Fermi liquid theory (FLT) [571] have also been employed in studies of low-energy charged- 
current neutrino-nucleus reactions. At finite temperature in stellar environments, thermal 
population of the excited states may enhance the weak interaction rates and cross sections 



at low neutrino energies 



58 



59|. 



Despite considerable progress in the development of advanced theoretical frameworks, up 
to this time only a limited number of microscopic models have been employed in large-scale 
calculations of neutrino-induced reactions and their implementation in supernova simula- 
tions. Because the complete modeling of neutrino-induced reactions necessitates the inclu- 
sion not only of Gamow- Teller transitions but also contributions from forbidden transitions 
and other higher multipoles, covering a large pool of target nuclei represents a computation- 
ally very demanding problem. Although the shell model provides a very accurate description 
of ground state wave functions, the description of high-lying excitations necessitates the use 
of large model spaces which often leads to computational difficulties, making the approach 
applicable essentially only to allowed transitions in light and medium-mass nuclei. For sys- 
tematic studies of neutrino-nucleus cross sections throughout the nuclide chart including the 
heavy nuclei, microscopic calculations must therefore be performed using models based on 
the QRPA. The first global calculations of electron neutrino-nucleus cross sections have been 



conducted using the gross theory of /3-decay 



60 



61] . Two microscopic frameworks have been 



employed in large-scale calculations of neutrino-induced reactions i) extended Thomas-Fermi 
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plus Strutinsky integral (ETFSI) and continuum quasiparticle random phase approximation 
(CQRPA) 62|, and ii) RPA with Landau-Migdal (LM) force, using Wood-Saxon potential 
(WS) to determine the single-particle basis of target nuclei [63j. There are no up-to-date 
systematic calculations conducted in the framework of energy density functional. 

Recently the framework based on relativistic nuclear energy density functional (RNEDF) 
has been introduced in modeling the charged-current neutrino- nucleus cross sections |64| . 
and it has also been extended for applications in modeling neutral-current neutrino-induced 
reactions 65|. In the case of iron group nuclei, comparison of the results for charge exchange 



reactions obtained using Skyrme functionals and the RNEDF, as well as with the shell model, 
showed reasonable theoretical uncertainty inherent in modeling neutrino-nucleus cross sec- 
tions. However, in view of applications in supernova and r-process simulations, which also 
involve unstable nuclei far from the valley of stability, where no experimental data are avail- 
able, it is necessary to provide independent insight into relevant neutrino-induced processes 
from various models and effective interactions. The main objective of this work is conduct- 
ing large-scale calculations of charged-current neutrino-nucleus cross sections in a large pool 
of nuclei from oxygen towards lead, based on the RNEDF. In addition to the overall cross 
sections covering the range of neutrino energies up to 100 MeV, calculations also include the 
cross sections averaged over neutrino fluxes for the range of temperatures characteristic of 
various stages of stellar evolution. In addition to the Fermi-Dirac neutrino spectra we also 
apply supernova spectra from a recent simulation that includes three-flavor Boltzmann neu- 
trino transport. Particular aim of this work is to emphasize the role of forbidden transitions 
in modeling neutrino-nucleus cross sections in large pool of nuclei, including both stable and 



nuclei away from the valley of stability. The re 



been discussed in several studies, e.g. Refs. 



24 



evance of forbidden transitions has already 



46 



66|. In the study of neutrino capture by 



r-process waiting point nuclei, first forbidden strength, together with the low-lying Gamow- 
Teller transitions, increased the rate of neutrino scattering from very neutron-rich nuclei by 
a factor of at least 2 and in some instances even by a factor of 5 |66|. As pointed out in mi- 
croscopic calculations based on Skyrme functionals in Ref. 46|, the properties of forbidden 
states are closely related to the neutrino-nucleus cross sections, and could be extracted by 
using neutrinos from low-energy beta-beams. 

The paper is organized as follows. Section II contains the basic theoretical background 
for the neutrino-nucleus cross sections in the charged-current channel based on weak Hamil- 
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tonian and the RNEDF. The results of large-scale calculations of the neutrino-nucleus cross 
sections are illustrated and discussed in Sec. III. The conclusions of the present work are 
summarized in Sec. IV. 



II. THEORETICAL BACKGROUND 



In the present work we explore the charged- current neutrino-nucleus process, 



(1) 



where the incoming electron neutrino (z/g) induces charge-exchange reaction in target nucleus 
X{Z,N). The general formalism of the neutrino-nucleus cross sections, derived assuming 



the Hami 
Refs. 53, 



tonian of the weak interaction in the current-current form, is given in details in 



54l |. The cross sections include the transition matrix elements between the initial 



I Jj) and final nuclear state \ Jf), for the charge Aij, longitudinal Cj, transverse electric Tf^, 
and transverse magnetic 7/'^"^*^ multipole operators 5J|. In the present work the RNEDF 
is employed in calculations of the transition matrix elements contributing to the neutrino- 
nucleus cross sections. The RNEDF has already been successfully employed in studies of 



giant resonances and exotic modes of excitation 



skin in nuclei 177 



68 



73], /3-decay rates of r-process nuclei 



muon capture l75| and stellar electron capture rates 76|, and in constraining the neutron 



More details about the implementation of the RNEDF in modeling 



64 



671. 



charged-current neutrino-nucleus reactions are given in Refs. 

The RNEDF based framework employs the self-consistent mean field for nucleons and 
minimal set of meson fields; isoscalar scalar cr-meson (J'^ = 0"*",T = 0), isoscalar vector u- 
meson {J^ = 1~,T = 0) and the isovector vector p-meson {J^ = 1~,T = 1), supplemented 
with the electromagnetic field. The meson-nucleon interaction is included with a minimal 
set of the interaction terms, where the vertex functionals include explicit dependence on 
the vector density. The nuclear ground state properties are described using the relativistic 
Hartree-Bogoliubov model (RHB), and relevant transitions induced by neutrinos are calcu- 
lated in the relativistic quasiparticle random phase approximation (RQRPA). More detai 
on the RHB model based on effective density-dependent interactions are given in Ref. 
The RQRPA is formulated in the canonical single-nucleon basis of the RHB model 



80 



79|. 



In modeling the neutrino-nucleus cross sections, important advantage of this framework is 
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that it is fully consistent in view of the effective interactions employed. In the particle-hole 
(ph) and pairing (pp) channels, the same interactions are used in the RHB equations that 
determine the canonical quasiparticle basis, and in the matrix equations of the RQRPA. In 
this way, one can employ the same RNEDF in description of the weak processes throughout 
the nuclide map without any additional adjustments to the properties of specific target nu- 
clei under consideration. For the model parameters that determine the density-dependent 



couplings and the meson masses, in this work DD-ME2 parameterization is used 



82|. The 



p^Hng correlafous in open shell nuclei are described by the finite range Gogrry interacfon. 

with parameterization DIS [83j . 

Complete calculation of inelastic neutrino-nucleus reactions spanning the range of neu- 
trino energies up to ^ 100 MeV necessitates the inclusion of a number of transitions with 
various multipoles J Although higher-order multipoles have rather small contributions 
at low incoming neutrino energies, these can not be neglected at energies about tens of 
MeV. In the present study, multipoles up to J = 5 have been taken into account. The large 
scale calculations of neutrino-nucleus cross sections, involving more than 1000 nuclei and a 
complete set of all multipoles up to J = 5 with both parities, necessitate considerable com- 
putational effort. Therefore, for the purposes of the present work, computational framework 
has been developed using parallel computing methods based on Message Passing Interface 



(MPI) SJl for the implementation on cluster and grid computer systems. 



III. RESULTS AND DISCUSSION 

By employing the model outlined in the previous section, we have conducted large-scale 
calculations of (z/g, e~) reactions in the OPb pool of even-even target nuclei spanning the 
range from oxygen toward lead (Z=8-82) with neutron number N=8-182. The calculations 
include excitations of all multipoles up to J = 5 and both parities. In the following we 
explore the systematic behavior of the overall cross sections throughout the pool under 
consideration, including the cross sections averaged over the Michel neutrino flux obtained 



from the decay at rest (DAR) of /x+ |l3| 



fmiE.) = -^{m,- . (2) 
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To simulate the supernova neutrino spectrum, at first instance we use the Fermi-Dirac 
distribution 

1 E"^ 

hii{Ey) = 7—— . (3) 

T-^ exp [[Ey/T) - 77J + 1 

In Fig. [1] the inclusive neutrino-nucleus cross sections, averaged over the Michel neutrino 
flux, are shown for the OPb pool of nuclei. Stable nuclei in the pool are denoted by flUed 
circles. In addition, the calculated cross sections for ^^C and ^^Fe are especially emphasized 



in comparison to the KARMEN experimental data |12|. One can observe that the pool 
of nuclei calculated with the same energy density functional without any adjustments of the 
model parameters flts nicely into two experimental data points. The cross sections increase 
systematically with increasing neutron number, resulting in values larger up to a factor 
^ 2 — 3 in comparison to those in the valley of stability. We note that sharp edge of the 
data at large A denotes boundary values for the pool under consideration, obtained for Pb 
isotope chain. 

Figure [2] shows the (z/e,e^) cross sections in the OPb pool of nuclei, averaged over the 
supernova neutrino flux given by Eq. ([3]) in the case of T = 4 MeV and 77 = 0. The results 
can be compared to those of the stable target nuclei. While the cross sections display a 
similar pattern as in the previous case with the Michel spectrum, the details of the data, 
however, depend on (T, 77) of the neutrino spectrum. The cross sections for various groups of 
target nuclei are separately displayed in Fig. [2) stable nuclei, proton-rich nuclei constrained 
by N/Z < 1, and neutron-rich nuclei with with N/Z > 1.5. The cross sections exhibit a 
systematic behavior throughout the nuclide map. In the case of neutron rich nuclei, the 
cross sections are larger in comparison to the stable nuclei, due to the increased number 
of neutrons that participate in charge-exchange neutrino-induced reactions. For proton-rich 
nuclei, the reaction pattern is opposite due to \)the smaller number of neutron-proton 2qp 
configurations, i.e., the {v^,e~) cross sections are considerably reduced in comparison to 
those of stable nuclei. This point is further illustrated in Fig. [3] where we plot the same 
cross sections versus the difference between the neutron and the proton numbers. For nuclei 
with N — Z > 1 the cross sections increase significantly with increasing number of excess 
neutrons. Although some scattering is apparent in the data, most of the results cluster 
along an almost linear function of the number of excess neutrons. Temperature dependence 
of neutrino-induced reactions within the OPb pool is illustrated in Fig. |U where the cross 
sections averaged over the Fermi-Dirac distribution (77 = 0) are shown for stable target nuclei 



in the range T = 2 — 8 MeV. At low temperatures, the cross sections appear rather scattered 
due to the stronger dependence on low-energy excitations in nuclei. However, for higher 
T, the cross sections are rather smooth and considerably larger due to the inclusion of a 
number of multipoles contributing at higher neutrino energies. The full set of calculations 
has been completed for the OPb pool in the range of temperatures T = — 10 MeV and 
complete tables are available on request. 

In view of applications in astrophysical models and in predicting the detector response to 
neutrinos involving various target nuclei, it is crucial to assess the sensitivity of the neutrino- 
nucleus cross sections on the theoretical frameworks and effective interactions employed. In 
the case of iron group nuclei, it has been shown that by employing different microscopic 
models, one can estimate reasonable theoretical uncertainty in neutrino-nucleus cross sec- 
tions averaged over the Michel spectrum, i.e., for ^'^Fe < a >th= (258 ±57) x IQ-'^^ cm^ j67|. 
Since for nuclei far from the valley of stability various models result in larger differences 
in nuclear structure properties, one could also expect larger sensitivity in the results of 
modeling neutrino-induced reactions. In the present analysis theoretical uncertainties are 
assessed on the basis of the OPb pool of nuclei. Figures |5l |6] and [7] show the dependence 
of the fiux-averaged neutrino-nucleus cross sections per nucleon (T = 4 MeV, r/ = 0) on the 
neutron number in the cases of Ni, Sn, and Pb isotopic chains, respectively. In these figures 
we compare our results to the available sets of cross sections based on two other theoretical 
frameworks: ETFSI + CQRPA |62| and RPA (WS+LM) joSl, where the former contains 
only contributions from the Gamow- Teller and Fermi transitions, while the latter includes 
all transitions up to J = 3. The results of the present RQRPA analysis for < a^^, > /A are 
shown separately for the full calculation including all transitions up to J = 5, and partial 
cross sections obtained only for the isobaric analog and Gamow- Teller transitions. 

The RQRPA results for N > Z nuclei, obtained taking into account only the Fermi 
and the Gamow- Teller transitions, are consistently higher than the corresponding values 
calculated using the ETFSI + CQRPA framework for all three isotopic chains. The deviation 
is rather small in the region around the valley of stability in Ni and Sn isotopes, but increases 
with additional neutrons. In the case of the Pb isotopic chain, the difference is significant 
even for the lightest isotopes considered in this study. For nuclei with the proton-to-neutron 
ratio Z/N < 1, the ETFSI + CQRPA framework predicts a change of the general trend and 
an anomalous increase of the cross sections with a corresponding decrease of the number of 
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neutrons. The agrement between results that take into accout forbidden transitions, RQRPA 
and RPA(WS+LM), is much better in all three isotopic chains, although for very neutron- 
rich Ni and Sn isotopes the RQRPA predicts higher values of the fiux-averaged cross section. 
In the case of the Pb isotopic chain, the agreement between the RQRPA and RPA(WS+LM) 
cross sections is excellent for all the isotopes under consideration. There are several reasons 
for the shown deviations. Each model employs different effective interactions that result 
in variations of the excitation pattern contributing to the cross sections, and the present 
calculations are based on a fully self-consistent approach to the neutrino-induced reactions. 
The three models agree best for nuclei around the valley of stability, where most of the 
experimental data on nuclear structure properties are available. In the very neutron rich 
region the differences between various models increase and grow relatively large. However, 
in all three isotopic chains the agreement of the RQRPA results is much better with the 
values obtained using the RPA (WS + LM) model than with the ETFSI + CQRPA, partly 
indicating the importance of the forbidden transitions. 

In the case of cross sections averaged over the supernova neutrino flux for (T = 4, rj = 0) , 
the impact of the self-consistent and complete calculations including all relevant multipoles 
is explored in view of previous knowledge on neutrino induced reactions on a large-scale basis 
within the OPb pool. Figure [8] shows the ratio of the cross sections averaged over supernova 



62l | model. In order to identify 



neutrino flux of the present work and ETFSI + CQRPA 
nuclei with pronounced discrepancies, the following groups of nuclei are separately denoted 
in figure: stable nuclei, proton-rich nuclei defined by N/Z < 1 and neutron-rich nuclei with 
N/Z > 1.5. For target nuclei where the cross sections are available for both models, the 
present large-scale calculations result mainly in systematically larger cross sections, for most 
of nuclei up to the factor of 4, and for smaller number of nuclei up to the factor of 7. However, 
for a limited set of medium mass nuclei the cross sections of this work are smaller than the 
ETFSI+CQRPA ones, and these are mainly limited to neutron deficient nuclei in the range 
50 ^ A ^ 100. The reason is the anomalous increase of the cross sections in the ETFSI + 
CQRPA model, when moving along a particular isotope chain from N = Z nucleus toward 
proton rich nuclei (see Figs. [5] and [6]). This behavior has not been observed in the present 
analysis. Apart from this anomaly, one can observe systematic linear increase of the ratio 
of cross sections with the nuclear mass. As already discussed in the case of isotopic chains, 
in addition to the self-consistent implementation of the RNEDF, the model employed in 
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the present work includes a complete set of transition operators up to J = 5 multipoles 
which enhance the overall cross sections. As shown in Fig. [8] the linear trend is particularly 
apparent in the case of neutron-rich nuclei {N/Z > 1.5), where the absolute values of the 
cross sections are relatively large in comparison to other nuclei under consideration. 

To further explore the role of forbidden transitions in the OPb pool, in Fig. IHlwe plot the 
ratio of cross sections taking into account all multipoles up to and including J = 5 with cross 
sections that only include the contributions from the Fermi and Gamow- Teller transitions. 
A trend very similar to the one in the previous figure appears, which indicates that the 
forbidden transitions contribute a larger fraction of the total cross sections in nuclei with 
higher mass. In heavier nuclei, and neutron-rich nuclei in particular, the difference between 
the proton and neutron numbers grows large, to the extent that in the very heavy nuclei 
neutrons occupy a full shell more above the protons. However, even in cases where there is 
not enough neutrons to occupy a full shell, additional neutrons move closer in energy to the 
orbits of the next shell, and therefore make the forbidden transitions less suppressed. The 
present analysis provides quantitative estimates for the mass dependence of the underlying 
structure in the neutrino-nucleus cross sections: in the mass region below A < 50 the 
forbidden transitions contribute less than 10% of the total flux- averaged cross sections, but 
their contribution grows practically linear with mass and can provide up to 50% of the 
total cross section. Therefore, in order to provide a realistic description of neutrino-nucleus 
reactions, as well as in the corresponding astrophysical applications, models must take into 
account more than the simplest Fermi and Gamow- Teller terms. 

In Fig. [To] the ratio of the averaged cross sections of this work and RPA (WS + LM) 63|] is 
shown. Qualitative agreement between the two models is obtained and differences are within 
a factor of 2. The ratio exhibits a mild mass dependence, i.e. the model based on RQRPA 
provides larger cross sections mainly for nuclei up to A ^ 150, while for heavier systems 
the cross sections are up to 50% smaller when compared to the ones previously reported 
with the RPA (WS + LM) model. These discrepancies occur for particular isotopic chains, 
however, there is no global systematic over- or underestimation of the cross-sections with 
respect to the RPA (WS + LM) approach, indicating that the deviations mainly originate 
from the different nuclear ground state and the residual interaction. The lack of a linear 
trend for neutron-rich nuclei (in comparison to Figs. [9] and [8]), along with the fact that RPA 
(WS+LM) model includes all transitions up to and including J = 3 points to the importance 
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of the lowest forbidden transitions in the quantitative description of semileptonic processes 
in nuclei. The largest deviations appear for low mass proton-rich nuclei, where the (z/g, e~) 
reactions become suppressed and the absolute values of the cross sections are actually very 
low, making them very sensitive to the structural details of the ground state. In comparison 
to the RPA (WS + LM) model, the main advantage of the present work is implementation 
of fully self-consistent theory framework in description of nuclear sector in neutrino induced 
reactions. In this way predictions for the cross sections become feasible not only in the region 
of stable nuclei and those that are covered by the RPA (WS+LM) model , but also in the 
unknown regions of the nuclide map of relevance for the stellar processes and nucleosynthesis. 

In addition to the Fermi-Dirac spectrum Eq. ([3]), we also apply the distribution obtained 
from a core-collapse supernova simulation (sSl. It can be described by the a-fit as follows 

where the parameter a is determined from the mean neutrino energy and the root-mean- 
square of the neutrino energy as |86| 



These data are taken from the core-collapse supernova explosion simulations of Ref. js?!. The 
supernova model is based on general relativistic radiation hydrodynamics and three-flavor 
Boltzmann neutrino transport in spherical symmetry. It also employs a nuclear equation of 



state [88|J and a modern set of neutrino opacities (for details, see Table 1 in Ref. jSSj and the 
references therein). The simulations were launched from massive progenitor stars, the 8.8 Mq 
0-Ne-Mg-core and the more massive Fe-core progenitors of 10.8 and 18 M0. Note that while 
for the low-mass star neutrino driven explosions can be obtained in spherically symmetric 
models, neutrino heating/cooling was enhanced artificially in order to trigger the onset of 



explosion for the more massive stars (for details, see Ref. (87|). For the current discussions, 
we select the data from the 18 Mq explosion model. Nevertheless, all simulations exhibit 
the same intrinsic feature of the continuously decreasing neutrino luminosities and average 
energies (see, e.g.. Figure 14 in Ref. js?!) after the supernova explosion has been launched. 
This important aspect is related to the continuous emission of neutrinos of all flavors, which 
deleptonizes the central protoneutron star on timescales of tens of seconds. Note that at the 
onset of the supernova explosion, the protoneutron star is hot (10-40 MeV) and still lepton 
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rich, in which sense it differs from the final supernova remnant neutron star. In order to 
characterize the dynamical evolution of the neutrino spectra during the protoneutron star 
deleptonization, we select three different post bounce times after the explosion has been 
launched, tpi, = 1, 5, and 20 seconds (see Table |T]for the corresponding neutrino energies). 

Once the neutrinos have decoupled from matter at the sphere of last scattering, which 
depends on neutrino flavor and energy, they can still contribute to charge and neutral cur- 
rent processes involving heavy nuclei at large distance (~ 10^ — 10^ km) relevant for the 



90] . At conditions where these 



nucleosynthesis, e.g., the up process |89[ and the r process 
processes can occur, the zero-temperature approach employed in the current paper for the 
calculation of the cross sections is a valid approximation. The protoneutron star deleptoniza- 
tion with the continuously reducing average neutrino energies has important consequences 
for the neutrino fiux-averaged cross sections, as shown in Fig. [TT] In comparison to the 
Fermi-Dirac spectrum with T = 4 MeV shown in Fig. |2l the supernova simulation spectrum 
results in significantly lower fiux-averaged cross sections by more than a factor of two for all 
nuclei under investigation, even shortly (tpb = 1 s) after the onset of explosion. Note that 
the supernova neutrino temperature, which can be defined as the matter temperature at the 
neutrinosphere, is slightly higher with T ^ 4.5 — 5 MeV than the assumed T = 4 MeV of the 
Fermi-Dirac spectrum. However, the Fermi-Dirac neutrino spectrum is significantly broader 
than those taken from the supernova simulations, which can be seen by comparing mean 
neutrino energies (E) and the mean-square energies {E'^) of /fd and /„. For comparison, 
we also plot the spectra in Fig. [121 The above mentioned difference between Fermi-Dirac 
spectra and supernova simulation spectra even increases during the ongoing explosion at late 
times, i.e. the supernova spectra shift increasingly towards Maxwell-Boltzmann like spectra. 
These findings are summarized in Tabled comparing the Fermi-Dirac spectra (/fd) and the 
supernova simulation spectra characterized by the a-fit (fa) at different post-bounce times 
after the explosion has been launched. 



IV. CONCLUSION 



We have employed a fully self-consistent model based on the RNEDF in large-scale calcu- 
lations of charged-current neutrino-nucleus cross sections in the OPb pool of target nuclei, 
spanning the range 8 < Z < 82 and 8 < < 182. The two main advantages of this ap- 
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TABLE I: Neutrino mean and mean-square energies for the different neutrino distributions. 



tpfe" [s] {E) [MeV] {E^) [MeV2] 



/fd 




12.60 


207.03 


fa 


1 


9.32 


108.57 


fa 


5 


8.68 


95.29 


fa 


20 


6.73 


58.32 



"Post bounce simulation time for the supernova spectra 

proach are: (i) self-consistent modeling of all relevant transition matrix elements involving 
open-shell nuclei, without any additional adjustments of the model parameters to the nu- 
clear target under consideration, and (ii) treating the allowed and the forbidden transitions 
on an equal footing, i.e., transitions of all relevant multipoles are taken into account. For the 
purpose of the present work, advanced computational framework has been developed, based 
on a parallel computing scheme using the Message Passing Interface for the implementation 
on cluster and grid computer systems. Reasonable agreement of the cross sections with 
the only available experimental data, for ^^C and ^^Fe, support the course of systematic 
calculations throughout the nuclide map. The results include a complete set of the cross 
sections for the OPb pool, calculated for the range of neutrino energies E^^ = — 100 MeV. 
The cross sections were averaged over the experimental Michel spectrum, the neutrino spec- 
trum described by the Fermi-Dirac distribution for the range of temperatures < T < 10 
MeV and chemical potential t] = 0. In addition, the averaged cross sections have been cal- 
culated using a-fit neutrino fluxes, taken from a recent core-collapse supernova explosion 
model that is based on three-flavor Boltzmann neutrino transport. Comparing Fermi-Dirac 
and supernova simulation neutrino spectra, we found that the latter result in significantly 
smaller fiux-averaged cross sections due to the more pinched supernova neutrino spectra. 
This aspect even increases during the ongoing deleptonization on the order of 10-20 s after 
which the cross sections become negligible and charged-current processes with nuclei have 
ceased. 

Note that the neutrino spectra form the supernova simulations do not take into account 
possible collective neutrino-fiavor oscillations. These phenomena can take place after neu- 
trino decoupling from matter and result in complete spectral swaps above a certain neutrino 
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energy for inverted neutrino-mass hierarchy 



9l| . Because heavy- lepton flavor neutrinos 



have higher average energies, being less strongly bound to matter in the absence of charged- 
current weak processes, it enhances the high-energy tail of the electron-flavor neutrino spec- 
tra. This jn turn may impact neutrino-induced nucleosynthesis of heavy elements at large 
distance 
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93|. 



Furthermore, it has been realized that medium modiflcations for the charged- current weak 
processes with neutrons and protons must be taken into account when computing neutrino 
transport and neutrino decoupling from matter at high densities using a nuclear equation 



of state [94|- These modiflcations have been explored at the mean-fleld level and shown to 



increase spectral differences between electron neutrinos and antineutrinos [95N97l|. Improved 
supernova simulations that take these effects into account may ultimately alter the results 
discussed in the current paper. 

Several key features of the neutrino-nucleus cross sections have been illustrated for various 
supernova neutrino fluxes, indicating systematic increase of the averaged cross sections with 
increase of the number of neutrons in target nuclei as well as with temperature of the 
neutrino distribution. When going away from the valley of stability toward neutron-rich 
nuclei, the cross sections become considerably enhanced, while on the proton-rich side they 
are suppressed due to the blocking of orbitals available for neutrino induced transitions. 

The cross sections from the present analysis have been discussed in view of previously 
reported large scale calculations. Current results are consistently higher than the results 
obtained using the ETFSI + CQRPA model which includes only the Fermi and the Gamow- 
Teller transitions. For lighter nuclei the RQRPA cross sections are approximately 50% 
larger, while for the heaviest nuclei studied they are up to 3 times larger. This enhancement 
originates in part from the forbidden transitions that are fully taken into account in the 
present work. We provide important quantitative insight into the underlying structure of 
the neutrino-nucleus cross sections throughout the nuclide map: in the mass region below 
A < 50 the forbidden transitions contribute less than 10% of the total flux-averaged cross 
sections, but their contribution increases with mass and can provide up to 50% of the total 
cross section. For the limited set of nuclei, the differences between the present results and 
those of the RPA (LM + WS) model are constrained within a factor of three. There is no 
general trend in the cross sections ratio between the two models, and the differences arise 
mainly due to the different approach in the calculations of the nuclear ground state and the 
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residual interaction in (Q)RPA. 

We have shown that the self-consistent calculations based on the RNEDF, which take 
into account transitions of all relevant multipoles, result in differences in comparison to 
previously reported cross sections, and for a considerable number of target nuclei the cross 
sections are larger. Therefore, revised calculations in modeling stellar evolution and r-process 
nucleosynthesis based on self-consistent descriptions of neutrino-induced reactions may al- 
low an updated insight into the origin of elements in the Universe. On the other hand, 
by employing various theory frameworks for neutrino-nucleus cross sections, one could esti- 
mate the uncertainties in the calculated element abundance patterns. The tables with the 
cross sections vs. neutrino energies of this work are available from the authors on request. 
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FIG. 1: (Color online) Inclusive neutrino- nucleus cross sections for target nuclei in the OPB pool, 
averaged over the Michel neutrino flux from muon DAR. The cross sections for stable nuclei are 
denoted by filled circles. The ernphasized results for ^^C and ^^Fe are shown in comparison to the 
experimental data (KARMEN) 
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FIG. 2: (Color online) Inclusive neutrino- nucleus cross sections for target nuclei in the OPb pool 
averaged over the Fermi-Dirac distribution for T = 4 MeV, r/ = 0. The cross sections for various 
groups of target nuclei are emphasized: stable nuclei (filled red circles), neutron-deficient nuclei 
with N jZ < 1 (filled black circles), neutron-rich nuclei with N/Z > 1.5 (solid squares). 
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FIG. 3: Inclusive neutrino-nucleus cross sections averaged over the Fermi-Dirac distribution with 
T = 4 MeV is plotted versus the difference between the neutron and the proton number. 
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FIG. 5: (Color online) Flux averaged supernova i/e-nucleus cross sections for Ni isotopes (T = 4 



MeV, rj = 0). The RHB+RQRPA cross sections calculated using al 
to those of ETFSI+CQRPA which include lAS+GT transitions 
cover multipoles up to J=3 [63 1. 



J < 5 transitions are compared 



631 and RPA (WS+LM) which 
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FIG. 6: (Color online) The same as in Fig. [5l but for Sn isotopic chain. 
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FIG. 7: (Color online) The same as in Fig. [5l but for Pb isotopic chain. 
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FIG. 8: (Color online) The ratio of the cross sections averaged over supernova neutrino flux (T = 4 



MeV, r/ = 0) calculated using the RHB+RQRPA and ETFSI+CQRPA [621 models. Various groups 
of nuclei are separately denoted: stable nuclei (solid circles), neutron-deficient nuclei with N jZ < 1 
(crosses), neutron-rich nuclei with N/Z > 1.5 (solid squares). 
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FIG. 9: (Color online) The ratio of the flux-averaged cross sections obtained by taking into account 
all multipole contributions and cross sections containing only Fermi and Gamow- Teller transitions. 
The Fermi-Dirac distribution was used for the neutrino spectrum with T = 4 MeV and r] = 0. 
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FIG. 10: (Color online) The ratio of the cross sections averaged over supernova neutrino distribution 
(T = 4 MeV, ?7 = 0) calculated using the RHB+RQRPA and RPA (WS+LM) E] models. Various 
groups of nuclei are separately denoted: stable nuclei (solid circles), neuton-deficient nuclei with 
N jZ < 1 (crosses), neutron-rich nuclei with N jZ > 1.5 (solid squares). 
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FIG. 11: (Color online) Flux-averaged cross sections of all nuclei in the OPb pool using the super- 
nova simulation spectra, comparing the different post-bounce times = 1, 5, 20 s. 




FIG. 12: (Color online) Fermi-Dirac /fd and supernova neutrino spectra /„. The latter correspond 
to the selected post-bounce times tpf, = 1, 5, 20 s. 
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